%
% Status : main Dynare file
%
% Warning : this file is generated automatically by Dynare
%           from model file (.mod)

if isoctave || matlab_ver_less_than('8.6')
    clear all
else
    clearvars -global
    clear_persistent_variables(fileparts(which('dynare')), false)
end
tic0 = tic;
% Define global variables.
global M_ options_ oo_ estim_params_ bayestopt_ dataset_ dataset_info estimation_info ys0_ ex0_
options_ = [];
M_.fname = 'model_file';
M_.dynare_version = '4.6.0';
oo_.dynare_version = '4.6.0';
options_.dynare_version = '4.6.0';
%
% Some global variables initialization
%
global_initialization;
diary off;
diary('model_file.log');
M_.exo_names = cell(6,1);
M_.exo_names_tex = cell(6,1);
M_.exo_names_long = cell(6,1);
M_.exo_names(1) = {'en'};
M_.exo_names_tex(1) = {'en'};
M_.exo_names_long(1) = {'en'};
M_.exo_names(2) = {'em'};
M_.exo_names_tex(2) = {'em'};
M_.exo_names_long(2) = {'em'};
M_.exo_names(3) = {'ee'};
M_.exo_names_tex(3) = {'ee'};
M_.exo_names_long(3) = {'ee'};
M_.exo_names(4) = {'eeta'};
M_.exo_names_tex(4) = {'eeta'};
M_.exo_names_long(4) = {'eeta'};
M_.exo_names(5) = {'ed'};
M_.exo_names_tex(5) = {'ed'};
M_.exo_names_long(5) = {'ed'};
M_.exo_names(6) = {'esigma'};
M_.exo_names_tex(6) = {'esigma'};
M_.exo_names_long(6) = {'esigma'};
M_.endo_names = cell(97,1);
M_.endo_names_tex = cell(97,1);
M_.endo_names_long = cell(97,1);
M_.endo_names(1) = {'uc'};
M_.endo_names_tex(1) = {'uc'};
M_.endo_names_long(1) = {'uc'};
M_.endo_names(2) = {'Euc'};
M_.endo_names_tex(2) = {'Euc'};
M_.endo_names_long(2) = {'Euc'};
M_.endo_names(3) = {'m'};
M_.endo_names_tex(3) = {'m'};
M_.endo_names_long(3) = {'m'};
M_.endo_names(4) = {'Em'};
M_.endo_names_tex(4) = {'Em'};
M_.endo_names_long(4) = {'Em'};
M_.endo_names(5) = {'mh'};
M_.endo_names_tex(5) = {'mh'};
M_.endo_names_long(5) = {'mh'};
M_.endo_names(6) = {'y'};
M_.endo_names_tex(6) = {'y'};
M_.endo_names_long(6) = {'y'};
M_.endo_names(7) = {'c'};
M_.endo_names_tex(7) = {'c'};
M_.endo_names_long(7) = {'c'};
M_.endo_names(8) = {'cc'};
M_.endo_names_tex(8) = {'cc'};
M_.endo_names_long(8) = {'cc'};
M_.endo_names(9) = {'cl'};
M_.endo_names_tex(9) = {'cl'};
M_.endo_names_long(9) = {'cl'};
M_.endo_names(10) = {'i'};
M_.endo_names_tex(10) = {'i'};
M_.endo_names_long(10) = {'i'};
M_.endo_names(11) = {'k'};
M_.endo_names_tex(11) = {'k'};
M_.endo_names_long(11) = {'k'};
M_.endo_names(12) = {'l'};
M_.endo_names_tex(12) = {'l'};
M_.endo_names_long(12) = {'l'};
M_.endo_names(13) = {'rk'};
M_.endo_names_tex(13) = {'rk'};
M_.endo_names_long(13) = {'rk'};
M_.endo_names(14) = {'w'};
M_.endo_names_tex(14) = {'w'};
M_.endo_names_long(14) = {'w'};
M_.endo_names(15) = {'pi'};
M_.endo_names_tex(15) = {'pi'};
M_.endo_names_long(15) = {'pi'};
M_.endo_names(16) = {'yn1'};
M_.endo_names_tex(16) = {'yn1'};
M_.endo_names_long(16) = {'yn1'};
M_.endo_names(17) = {'ph'};
M_.endo_names_tex(17) = {'ph'};
M_.endo_names_long(17) = {'ph'};
M_.endo_names(18) = {'pe'};
M_.endo_names_tex(18) = {'pe'};
M_.endo_names_long(18) = {'pe'};
M_.endo_names(19) = {'pih'};
M_.endo_names_tex(19) = {'pih'};
M_.endo_names_long(19) = {'pih'};
M_.endo_names(20) = {'pie'};
M_.endo_names_tex(20) = {'pie'};
M_.endo_names_long(20) = {'pie'};
M_.endo_names(21) = {'pihnum'};
M_.endo_names_tex(21) = {'pihnum'};
M_.endo_names_long(21) = {'pihnum'};
M_.endo_names(22) = {'lambda'};
M_.endo_names_tex(22) = {'lambda'};
M_.endo_names_long(22) = {'lambda'};
M_.endo_names(23) = {'qk'};
M_.endo_names_tex(23) = {'qk'};
M_.endo_names_long(23) = {'qk'};
M_.endo_names(24) = {'a'};
M_.endo_names_tex(24) = {'a'};
M_.endo_names_long(24) = {'a'};
M_.endo_names(25) = {'dn'};
M_.endo_names_tex(25) = {'dn'};
M_.endo_names_long(25) = {'dn'};
M_.endo_names(26) = {'sigmaa'};
M_.endo_names_tex(26) = {'sigmaa'};
M_.endo_names_long(26) = {'sigmaa'};
M_.endo_names(27) = {'eta'};
M_.endo_names_tex(27) = {'eta'};
M_.endo_names_long(27) = {'eta'};
M_.endo_names(28) = {'ce'};
M_.endo_names_tex(28) = {'ce'};
M_.endo_names_long(28) = {'ce'};
M_.endo_names(29) = {'epe'};
M_.endo_names_tex(29) = {'epe'};
M_.endo_names_long(29) = {'epe'};
M_.endo_names(30) = {'etalag'};
M_.endo_names_tex(30) = {'etalag'};
M_.endo_names_long(30) = {'etalag'};
M_.endo_names(31) = {'dc'};
M_.endo_names_tex(31) = {'dc'};
M_.endo_names_long(31) = {'dc'};
M_.endo_names(32) = {'dcl'};
M_.endo_names_tex(32) = {'dcl'};
M_.endo_names_long(32) = {'dcl'};
M_.endo_names(33) = {'dcc'};
M_.endo_names_tex(33) = {'dcc'};
M_.endo_names_long(33) = {'dcc'};
M_.endo_names(34) = {'dce'};
M_.endo_names_tex(34) = {'dce'};
M_.endo_names_long(34) = {'dce'};
M_.endo_names(35) = {'dcr'};
M_.endo_names_tex(35) = {'dcr'};
M_.endo_names_long(35) = {'dcr'};
M_.endo_names(36) = {'dz'};
M_.endo_names_tex(36) = {'dz'};
M_.endo_names_long(36) = {'dz'};
M_.endo_names(37) = {'Edz'};
M_.endo_names_tex(37) = {'Edz'};
M_.endo_names_long(37) = {'Edz'};
M_.endo_names(38) = {'dy'};
M_.endo_names_tex(38) = {'dy'};
M_.endo_names_long(38) = {'dy'};
M_.endo_names(39) = {'di'};
M_.endo_names_tex(39) = {'di'};
M_.endo_names_long(39) = {'di'};
M_.endo_names(40) = {'Phik'};
M_.endo_names_tex(40) = {'Phik'};
M_.endo_names_long(40) = {'Phik'};
M_.endo_names(41) = {'Phik1'};
M_.endo_names_tex(41) = {'Phik1'};
M_.endo_names_long(41) = {'Phik1'};
M_.endo_names(42) = {'Fk'};
M_.endo_names_tex(42) = {'Fk'};
M_.endo_names_long(42) = {'Fk'};
M_.endo_names(43) = {'Fl'};
M_.endo_names_tex(43) = {'Fl'};
M_.endo_names_long(43) = {'Fl'};
M_.endo_names(44) = {'Jk'};
M_.endo_names_tex(44) = {'Jk'};
M_.endo_names_long(44) = {'Jk'};
M_.endo_names(45) = {'Jl'};
M_.endo_names_tex(45) = {'Jl'};
M_.endo_names_long(45) = {'Jl'};
M_.endo_names(46) = {'Phip'};
M_.endo_names_tex(46) = {'Phip'};
M_.endo_names_long(46) = {'Phip'};
M_.endo_names(47) = {'Phip1'};
M_.endo_names_tex(47) = {'Phip1'};
M_.endo_names_long(47) = {'Phip1'};
M_.endo_names(48) = {'Phip2'};
M_.endo_names_tex(48) = {'Phip2'};
M_.endo_names_long(48) = {'Phip2'};
M_.endo_names(49) = {'rf'};
M_.endo_names_tex(49) = {'rf'};
M_.endo_names_long(49) = {'rf'};
M_.endo_names(50) = {'d'};
M_.endo_names_tex(50) = {'d'};
M_.endo_names_long(50) = {'d'};
M_.endo_names(51) = {'dagg'};
M_.endo_names_tex(51) = {'dagg'};
M_.endo_names_long(51) = {'dagg'};
M_.endo_names(52) = {'dd'};
M_.endo_names_tex(52) = {'dd'};
M_.endo_names_long(52) = {'dd'};
M_.endo_names(53) = {'ddn'};
M_.endo_names_tex(53) = {'ddn'};
M_.endo_names_long(53) = {'ddn'};
M_.endo_names(54) = {'ddr'};
M_.endo_names_tex(54) = {'ddr'};
M_.endo_names_long(54) = {'ddr'};
M_.endo_names(55) = {'ddagg'};
M_.endo_names_tex(55) = {'ddagg'};
M_.endo_names_long(55) = {'ddagg'};
M_.endo_names(56) = {'ddaggn'};
M_.endo_names_tex(56) = {'ddaggn'};
M_.endo_names_long(56) = {'ddaggn'};
M_.endo_names(57) = {'ddaggr'};
M_.endo_names_tex(57) = {'ddaggr'};
M_.endo_names_long(57) = {'ddaggr'};
M_.endo_names(58) = {'pd'};
M_.endo_names_tex(58) = {'pd'};
M_.endo_names_long(58) = {'pd'};
M_.endo_names(59) = {'pdagg'};
M_.endo_names_tex(59) = {'pdagg'};
M_.endo_names_long(59) = {'pdagg'};
M_.endo_names(60) = {'rd'};
M_.endo_names_tex(60) = {'rd'};
M_.endo_names_long(60) = {'rd'};
M_.endo_names(61) = {'rdr'};
M_.endo_names_tex(61) = {'rdr'};
M_.endo_names_long(61) = {'rdr'};
M_.endo_names(62) = {'rdagg'};
M_.endo_names_tex(62) = {'rdagg'};
M_.endo_names_long(62) = {'rdagg'};
M_.endo_names(63) = {'rdaggn'};
M_.endo_names_tex(63) = {'rdaggn'};
M_.endo_names_long(63) = {'rdaggn'};
M_.endo_names(64) = {'rdaggr'};
M_.endo_names_tex(64) = {'rdaggr'};
M_.endo_names_long(64) = {'rdaggr'};
M_.endo_names(65) = {'rm'};
M_.endo_names_tex(65) = {'rm'};
M_.endo_names_long(65) = {'rm'};
M_.endo_names(66) = {'rmex'};
M_.endo_names_tex(66) = {'rmex'};
M_.endo_names_long(66) = {'rmex'};
M_.endo_names(67) = {'yn2'};
M_.endo_names_tex(67) = {'yn2'};
M_.endo_names_long(67) = {'yn2'};
M_.endo_names(68) = {'yn3'};
M_.endo_names_tex(68) = {'yn3'};
M_.endo_names_long(68) = {'yn3'};
M_.endo_names(69) = {'yn4'};
M_.endo_names_tex(69) = {'yn4'};
M_.endo_names_long(69) = {'yn4'};
M_.endo_names(70) = {'yn5'};
M_.endo_names_tex(70) = {'yn5'};
M_.endo_names_long(70) = {'yn5'};
M_.endo_names(71) = {'yn6'};
M_.endo_names_tex(71) = {'yn6'};
M_.endo_names_long(71) = {'yn6'};
M_.endo_names(72) = {'yn7'};
M_.endo_names_tex(72) = {'yn7'};
M_.endo_names_long(72) = {'yn7'};
M_.endo_names(73) = {'yn8'};
M_.endo_names_tex(73) = {'yn8'};
M_.endo_names_long(73) = {'yn8'};
M_.endo_names(74) = {'yn9'};
M_.endo_names_tex(74) = {'yn9'};
M_.endo_names_long(74) = {'yn9'};
M_.endo_names(75) = {'yn10'};
M_.endo_names_tex(75) = {'yn10'};
M_.endo_names_long(75) = {'yn10'};
M_.endo_names(76) = {'yn11'};
M_.endo_names_tex(76) = {'yn11'};
M_.endo_names_long(76) = {'yn11'};
M_.endo_names(77) = {'yn12'};
M_.endo_names_tex(77) = {'yn12'};
M_.endo_names_long(77) = {'yn12'};
M_.endo_names(78) = {'yn13'};
M_.endo_names_tex(78) = {'yn13'};
M_.endo_names_long(78) = {'yn13'};
M_.endo_names(79) = {'yn14'};
M_.endo_names_tex(79) = {'yn14'};
M_.endo_names_long(79) = {'yn14'};
M_.endo_names(80) = {'yn15'};
M_.endo_names_tex(80) = {'yn15'};
M_.endo_names_long(80) = {'yn15'};
M_.endo_names(81) = {'yn16'};
M_.endo_names_tex(81) = {'yn16'};
M_.endo_names_long(81) = {'yn16'};
M_.endo_names(82) = {'yn17'};
M_.endo_names_tex(82) = {'yn17'};
M_.endo_names_long(82) = {'yn17'};
M_.endo_names(83) = {'yn18'};
M_.endo_names_tex(83) = {'yn18'};
M_.endo_names_long(83) = {'yn18'};
M_.endo_names(84) = {'yn19'};
M_.endo_names_tex(84) = {'yn19'};
M_.endo_names_long(84) = {'yn19'};
M_.endo_names(85) = {'yn20'};
M_.endo_names_tex(85) = {'yn20'};
M_.endo_names_long(85) = {'yn20'};
M_.endo_names(86) = {'slope'};
M_.endo_names_tex(86) = {'slope'};
M_.endo_names_long(86) = {'slope'};
M_.endo_names(87) = {'fen'};
M_.endo_names_tex(87) = {'fen'};
M_.endo_names_long(87) = {'fen'};
M_.endo_names(88) = {'cf'};
M_.endo_names_tex(88) = {'cf'};
M_.endo_names_long(88) = {'cf'};
M_.endo_names(89) = {'dr'};
M_.endo_names_tex(89) = {'dr'};
M_.endo_names_long(89) = {'dr'};
M_.endo_names(90) = {'cfn'};
M_.endo_names_tex(90) = {'cfn'};
M_.endo_names_long(90) = {'cfn'};
M_.endo_names(91) = {'drn'};
M_.endo_names_tex(91) = {'drn'};
M_.endo_names_long(91) = {'drn'};
M_.endo_names(92) = {'drrf'};
M_.endo_names_tex(92) = {'drrf'};
M_.endo_names_long(92) = {'drrf'};
M_.endo_names(93) = {'dry'};
M_.endo_names_tex(93) = {'dry'};
M_.endo_names_long(93) = {'dry'};
M_.endo_names(94) = {'rmagg'};
M_.endo_names_tex(94) = {'rmagg'};
M_.endo_names_long(94) = {'rmagg'};
M_.endo_names(95) = {'rmexagg'};
M_.endo_names_tex(95) = {'rmexagg'};
M_.endo_names_long(95) = {'rmexagg'};
M_.endo_names(96) = {'Erd'};
M_.endo_names_tex(96) = {'Erd'};
M_.endo_names_long(96) = {'Erd'};
M_.endo_names(97) = {'Erdn'};
M_.endo_names_tex(97) = {'Erdn'};
M_.endo_names_long(97) = {'Erdn'};
M_.endo_partitions = struct();
M_.param_names = cell(64,1);
M_.param_names_tex = cell(64,1);
M_.param_names_long = cell(64,1);
M_.param_names(1) = {'delta'};
M_.param_names_tex(1) = {'delta'};
M_.param_names_long(1) = {'delta'};
M_.param_names(2) = {'psi'};
M_.param_names_tex(2) = {'psi'};
M_.param_names_long(2) = {'psi'};
M_.param_names(3) = {'gamma'};
M_.param_names_tex(3) = {'gamma'};
M_.param_names_long(3) = {'gamma'};
M_.param_names(4) = {'xil'};
M_.param_names_tex(4) = {'xil'};
M_.param_names_long(4) = {'xil'};
M_.param_names(5) = {'omega'};
M_.param_names_tex(5) = {'omega'};
M_.param_names_long(5) = {'omega'};
M_.param_names(6) = {'lbar'};
M_.param_names_tex(6) = {'lbar'};
M_.param_names_long(6) = {'lbar'};
M_.param_names(7) = {'l_ss'};
M_.param_names_tex(7) = {'l\_ss'};
M_.param_names_long(7) = {'l_ss'};
M_.param_names(8) = {'mu'};
M_.param_names_tex(8) = {'mu'};
M_.param_names_long(8) = {'mu'};
M_.param_names(9) = {'alphac'};
M_.param_names_tex(9) = {'alphac'};
M_.param_names_long(9) = {'alphac'};
M_.param_names(10) = {'phic'};
M_.param_names_tex(10) = {'phic'};
M_.param_names_long(10) = {'phic'};
M_.param_names(11) = {'nu'};
M_.param_names_tex(11) = {'nu'};
M_.param_names_long(11) = {'nu'};
M_.param_names(12) = {'alpha'};
M_.param_names_tex(12) = {'alpha'};
M_.param_names_long(12) = {'alpha'};
M_.param_names(13) = {'deltak'};
M_.param_names_tex(13) = {'deltak'};
M_.param_names_long(13) = {'deltak'};
M_.param_names(14) = {'phip'};
M_.param_names_tex(14) = {'phip'};
M_.param_names_long(14) = {'phip'};
M_.param_names(15) = {'xik'};
M_.param_names_tex(15) = {'xik'};
M_.param_names_long(15) = {'xik'};
M_.param_names(16) = {'rhon'};
M_.param_names_tex(16) = {'rhon'};
M_.param_names_long(16) = {'rhon'};
M_.param_names(17) = {'phin'};
M_.param_names_tex(17) = {'phin'};
M_.param_names_long(17) = {'phin'};
M_.param_names(18) = {'a1k'};
M_.param_names_tex(18) = {'a1k'};
M_.param_names_long(18) = {'a1k'};
M_.param_names(19) = {'a2k'};
M_.param_names_tex(19) = {'a2k'};
M_.param_names_long(19) = {'a2k'};
M_.param_names(20) = {'philev'};
M_.param_names_tex(20) = {'philev'};
M_.param_names_long(20) = {'philev'};
M_.param_names(21) = {'rhopd'};
M_.param_names_tex(21) = {'rhopd'};
M_.param_names_long(21) = {'rhopd'};
M_.param_names(22) = {'pd_ss'};
M_.param_names_tex(22) = {'pd\_ss'};
M_.param_names_long(22) = {'pd_ss'};
M_.param_names(23) = {'phia'};
M_.param_names_tex(23) = {'phia'};
M_.param_names_long(23) = {'phia'};
M_.param_names(24) = {'rhoa'};
M_.param_names_tex(24) = {'rhoa'};
M_.param_names_long(24) = {'rhoa'};
M_.param_names(25) = {'a0'};
M_.param_names_tex(25) = {'a0'};
M_.param_names_long(25) = {'a0'};
M_.param_names(26) = {'rhosigma'};
M_.param_names_tex(26) = {'rhosigma'};
M_.param_names_long(26) = {'rhosigma'};
M_.param_names(27) = {'sigmaa0'};
M_.param_names_tex(27) = {'sigmaa0'};
M_.param_names_long(27) = {'sigmaa0'};
M_.param_names(28) = {'sigmasigma'};
M_.param_names_tex(28) = {'sigmasigma'};
M_.param_names_long(28) = {'sigmasigma'};
M_.param_names(29) = {'sigmatarget'};
M_.param_names_tex(29) = {'sigmatarget'};
M_.param_names_long(29) = {'sigmatarget'};
M_.param_names(30) = {'rhosigma_a'};
M_.param_names_tex(30) = {'rhosigma\_a'};
M_.param_names_long(30) = {'rhosigma_a'};
M_.param_names(31) = {'corr_vol'};
M_.param_names_tex(31) = {'corr\_vol'};
M_.param_names_long(31) = {'corr_vol'};
M_.param_names(32) = {'sigmad'};
M_.param_names_tex(32) = {'sigmad'};
M_.param_names_long(32) = {'sigmad'};
M_.param_names(33) = {'sigmae'};
M_.param_names_tex(33) = {'sigmae'};
M_.param_names_long(33) = {'sigmae'};
M_.param_names(34) = {'sigmaeta'};
M_.param_names_tex(34) = {'sigmaeta'};
M_.param_names_long(34) = {'sigmaeta'};
M_.param_names(35) = {'rhoe'};
M_.param_names_tex(35) = {'rhoe'};
M_.param_names_long(35) = {'rhoe'};
M_.param_names(36) = {'rhoeta1'};
M_.param_names_tex(36) = {'rhoeta1'};
M_.param_names_long(36) = {'rhoeta1'};
M_.param_names(37) = {'rhoeta2'};
M_.param_names_tex(37) = {'rhoeta2'};
M_.param_names_long(37) = {'rhoeta2'};
M_.param_names(38) = {'deltaa'};
M_.param_names_tex(38) = {'deltaa'};
M_.param_names_long(38) = {'deltaa'};
M_.param_names(39) = {'ecap'};
M_.param_names_tex(39) = {'ecap'};
M_.param_names_long(39) = {'ecap'};
M_.param_names(40) = {'rhor'};
M_.param_names_tex(40) = {'rhor'};
M_.param_names_long(40) = {'rhor'};
M_.param_names(41) = {'rhopi'};
M_.param_names_tex(41) = {'rhopi'};
M_.param_names_long(41) = {'rhopi'};
M_.param_names(42) = {'rhoy'};
M_.param_names_tex(42) = {'rhoy'};
M_.param_names_long(42) = {'rhoy'};
M_.param_names(43) = {'sigmam'};
M_.param_names_tex(43) = {'sigmam'};
M_.param_names_long(43) = {'sigmam'};
M_.param_names(44) = {'pi_ss'};
M_.param_names_tex(44) = {'pi\_ss'};
M_.param_names_long(44) = {'pi_ss'};
M_.param_names(45) = {'yn1_ss'};
M_.param_names_tex(45) = {'yn1\_ss'};
M_.param_names_long(45) = {'yn1_ss'};
M_.param_names(46) = {'y_ss'};
M_.param_names_tex(46) = {'y\_ss'};
M_.param_names_long(46) = {'y_ss'};
M_.param_names(47) = {'cc_ss'};
M_.param_names_tex(47) = {'cc\_ss'};
M_.param_names_long(47) = {'cc_ss'};
M_.param_names(48) = {'ce_ss'};
M_.param_names_tex(48) = {'ce\_ss'};
M_.param_names_long(48) = {'ce_ss'};
M_.param_names(49) = {'c_ss'};
M_.param_names_tex(49) = {'c\_ss'};
M_.param_names_long(49) = {'c_ss'};
M_.param_names(50) = {'k_ss'};
M_.param_names_tex(50) = {'k\_ss'};
M_.param_names_long(50) = {'k_ss'};
M_.param_names(51) = {'i_ss'};
M_.param_names_tex(51) = {'i\_ss'};
M_.param_names_long(51) = {'i_ss'};
M_.param_names(52) = {'m_ss'};
M_.param_names_tex(52) = {'m\_ss'};
M_.param_names_long(52) = {'m_ss'};
M_.param_names(53) = {'r_ss'};
M_.param_names_tex(53) = {'r\_ss'};
M_.param_names_long(53) = {'r_ss'};
M_.param_names(54) = {'rf_ss'};
M_.param_names_tex(54) = {'rf\_ss'};
M_.param_names_long(54) = {'rf_ss'};
M_.param_names(55) = {'Fl_ss'};
M_.param_names_tex(55) = {'Fl\_ss'};
M_.param_names_long(55) = {'Fl_ss'};
M_.param_names(56) = {'Fk_ss'};
M_.param_names_tex(56) = {'Fk\_ss'};
M_.param_names_long(56) = {'Fk_ss'};
M_.param_names(57) = {'Jl_ss'};
M_.param_names_tex(57) = {'Jl\_ss'};
M_.param_names_long(57) = {'Jl_ss'};
M_.param_names(58) = {'Jk_ss'};
M_.param_names_tex(58) = {'Jk\_ss'};
M_.param_names_long(58) = {'Jk_ss'};
M_.param_names(59) = {'w_ss'};
M_.param_names_tex(59) = {'w\_ss'};
M_.param_names_long(59) = {'w_ss'};
M_.param_names(60) = {'qk_ss'};
M_.param_names_tex(60) = {'qk\_ss'};
M_.param_names_long(60) = {'qk_ss'};
M_.param_names(61) = {'pe_ss'};
M_.param_names_tex(61) = {'pe\_ss'};
M_.param_names_long(61) = {'pe_ss'};
M_.param_names(62) = {'Phik_ss'};
M_.param_names_tex(62) = {'Phik\_ss'};
M_.param_names_long(62) = {'Phik_ss'};
M_.param_names(63) = {'Phik1_ss'};
M_.param_names_tex(63) = {'Phik1\_ss'};
M_.param_names_long(63) = {'Phik1_ss'};
M_.param_names(64) = {'ph_ss'};
M_.param_names_tex(64) = {'ph\_ss'};
M_.param_names_long(64) = {'ph_ss'};
M_.param_partitions = struct();
M_.exo_det_nbr = 0;
M_.exo_nbr = 6;
M_.endo_nbr = 97;
M_.param_nbr = 64;
M_.orig_endo_nbr = 97;
M_.aux_vars = [];
M_.Sigma_e = zeros(6, 6);
M_.Correlation_matrix = eye(6, 6);
M_.H = 0;
M_.Correlation_matrix_ME = 1;
M_.sigma_e_is_diagonal = true;
M_.det_shocks = [];
options_.linear = false;
options_.block = false;
options_.bytecode = false;
options_.use_dll = false;
options_.linear_decomposition = false;
M_.nonzero_hessian_eqs = [1 2 3 4 5 6 7 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 39 40 41 44 47 48 49 50 51 52 53 54 55 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 80 81 82 85 86 88 89];
M_.hessian_eq_zero = isempty(M_.nonzero_hessian_eqs);
M_.orig_eq_nbr = 97;
M_.eq_nbr = 97;
M_.ramsey_eq_nbr = 0;
M_.set_auxiliary_variables = exist(['./+' M_.fname '/set_auxiliary_variables.m'], 'file') == 2;
M_.epilogue_names = {};
M_.epilogue_var_list_ = {};
M_.orig_maximum_endo_lag = 1;
M_.orig_maximum_endo_lead = 1;
M_.orig_maximum_exo_lag = 0;
M_.orig_maximum_exo_lead = 0;
M_.orig_maximum_exo_det_lag = 0;
M_.orig_maximum_exo_det_lead = 0;
M_.orig_maximum_lag = 1;
M_.orig_maximum_lead = 1;
M_.orig_maximum_lag_with_diffs_expanded = 1;
M_.lead_lag_incidence = [
 0 23 120;
 1 24 0;
 0 25 121;
 0 26 0;
 0 27 0;
 2 28 0;
 3 29 0;
 4 30 0;
 5 31 0;
 6 32 0;
 7 33 0;
 0 34 0;
 0 35 122;
 0 36 0;
 0 37 123;
 8 38 124;
 9 39 0;
 10 40 125;
 0 41 0;
 0 42 0;
 0 43 126;
 0 44 0;
 11 45 0;
 12 46 0;
 13 47 0;
 14 48 0;
 15 49 0;
 16 50 0;
 0 51 0;
 17 52 0;
 0 53 0;
 0 54 127;
 0 55 0;
 0 56 0;
 0 57 0;
 0 58 128;
 0 59 0;
 0 60 0;
 0 61 0;
 0 62 0;
 0 63 0;
 0 64 0;
 0 65 0;
 0 66 0;
 0 67 0;
 0 68 0;
 0 69 0;
 0 70 129;
 18 71 0;
 19 72 0;
 20 73 0;
 0 74 130;
 0 75 0;
 0 76 0;
 0 77 131;
 0 78 132;
 0 79 133;
 21 80 134;
 22 81 135;
 0 82 0;
 0 83 0;
 0 84 136;
 0 85 137;
 0 86 138;
 0 87 139;
 0 88 0;
 0 89 140;
 0 90 141;
 0 91 142;
 0 92 143;
 0 93 144;
 0 94 145;
 0 95 146;
 0 96 147;
 0 97 148;
 0 98 149;
 0 99 150;
 0 100 151;
 0 101 152;
 0 102 153;
 0 103 154;
 0 104 155;
 0 105 156;
 0 106 157;
 0 107 0;
 0 108 0;
 0 109 0;
 0 110 158;
 0 111 159;
 0 112 160;
 0 113 161;
 0 114 162;
 0 115 163;
 0 116 164;
 0 117 0;
 0 118 0;
 0 119 0;]';
M_.nstatic = 34;
M_.nfwrd   = 41;
M_.npred   = 18;
M_.nboth   = 4;
M_.nsfwrd   = 45;
M_.nspred   = 22;
M_.ndynamic   = 63;
M_.dynamic_tmp_nbr = [43; 56; 52; 16; ];
M_.equations_tags = {
  1 , 'name' , '1' ;
  2 , 'name' , '2' ;
  3 , 'name' , 'mh' ;
  4 , 'name' , '4' ;
  5 , 'name' , 'm' ;
  6 , 'name' , '6' ;
  7 , 'name' , 'y' ;
  8 , 'name' , '8' ;
  9 , 'name' , '9' ;
  10 , 'name' , 'Phik' ;
  11 , 'name' , '11' ;
  12 , 'name' , 'Fk' ;
  13 , 'name' , 'Jk' ;
  14 , 'name' , 'Phik1' ;
  15 , 'name' , '15' ;
  16 , 'name' , '16' ;
  17 , 'name' , '17' ;
  18 , 'name' , 'Fl' ;
  19 , 'name' , 'Jl' ;
  20 , 'name' , '20' ;
  21 , 'name' , 'Phip' ;
  22 , 'name' , 'Phip1' ;
  23 , 'name' , 'Phip2' ;
  24 , 'name' , '24' ;
  25 , 'name' , '25' ;
  26 , 'name' , '26' ;
  27 , 'name' , 'pe' ;
  28 , 'name' , '28' ;
  29 , 'name' , 'pie' ;
  30 , 'name' , 'pih' ;
  31 , 'name' , '31' ;
  32 , 'name' , 'epe' ;
  33 , 'name' , 'a' ;
  34 , 'name' , 'dn' ;
  35 , 'name' , '35' ;
  36 , 'name' , 'etalag' ;
  37 , 'name' , 'eta' ;
  38 , 'name' , 'sigmaa' ;
  39 , 'name' , 'd' ;
  40 , 'name' , 'dagg' ;
  41 , 'name' , 'dd' ;
  42 , 'name' , 'ddr' ;
  43 , 'name' , 'ddn' ;
  44 , 'name' , 'ddagg' ;
  45 , 'name' , 'ddaggr' ;
  46 , 'name' , 'ddaggn' ;
  47 , 'name' , '47' ;
  48 , 'name' , '48' ;
  49 , 'name' , '49' ;
  50 , 'name' , '50' ;
  51 , 'name' , '51' ;
  52 , 'name' , '52' ;
  53 , 'name' , '53' ;
  54 , 'name' , 'rm' ;
  55 , 'name' , 'rmagg' ;
  56 , 'name' , 'rmex' ;
  57 , 'name' , 'rmexagg' ;
  58 , 'name' , 'Em' ;
  59 , 'name' , '59' ;
  60 , 'name' , '60' ;
  61 , 'name' , '61' ;
  62 , 'name' , '62' ;
  63 , 'name' , '63' ;
  64 , 'name' , '64' ;
  65 , 'name' , '65' ;
  66 , 'name' , '66' ;
  67 , 'name' , '67' ;
  68 , 'name' , '68' ;
  69 , 'name' , '69' ;
  70 , 'name' , '70' ;
  71 , 'name' , '71' ;
  72 , 'name' , '72' ;
  73 , 'name' , '73' ;
  74 , 'name' , '74' ;
  75 , 'name' , '75' ;
  76 , 'name' , '76' ;
  77 , 'name' , '77' ;
  78 , 'name' , '78' ;
  79 , 'name' , 'slope' ;
  80 , 'name' , 'dc' ;
  81 , 'name' , 'dcc' ;
  82 , 'name' , 'dce' ;
  83 , 'name' , 'dz' ;
  84 , 'name' , 'Edz' ;
  85 , 'name' , 'cl' ;
  86 , 'name' , 'dcl' ;
  87 , 'name' , 'dcr' ;
  88 , 'name' , 'dy' ;
  89 , 'name' , 'di' ;
  90 , 'name' , 'cf' ;
  91 , 'name' , 'dr' ;
  92 , 'name' , 'cfn' ;
  93 , 'name' , 'drn' ;
  94 , 'name' , 'drrf' ;
  95 , 'name' , 'dry' ;
  96 , 'name' , 'Erd' ;
  97 , 'name' , 'Erdn' ;
};
M_.mapping.uc.eqidx = [1 2 3 ];
M_.mapping.Euc.eqidx = [1 2 3 ];
M_.mapping.m.eqidx = [5 15 20 24 26 47 50 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 ];
M_.mapping.Em.eqidx = [58 ];
M_.mapping.mh.eqidx = [3 5 ];
M_.mapping.y.eqidx = [7 8 12 18 21 22 23 25 39 88 ];
M_.mapping.c.eqidx = [4 6 28 80 85 ];
M_.mapping.cc.eqidx = [4 8 27 28 81 ];
M_.mapping.cl.eqidx = [85 86 ];
M_.mapping.i.eqidx = [8 10 11 14 39 89 ];
M_.mapping.k.eqidx = [7 9 10 11 12 13 14 ];
M_.mapping.l.eqidx = [6 7 18 19 39 85 ];
M_.mapping.rk.eqidx = [11 15 ];
M_.mapping.w.eqidx = [6 17 39 ];
M_.mapping.pi.eqidx = [21 22 23 24 25 26 29 30 31 42 45 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 88 89 ];
M_.mapping.yn1.eqidx = [24 25 60 79 95 ];
M_.mapping.ph.eqidx = [5 6 28 30 ];
M_.mapping.pe.eqidx = [27 28 29 31 32 40 59 88 ];
M_.mapping.pih.eqidx = [30 87 ];
M_.mapping.pie.eqidx = [29 ];
M_.mapping.pihnum.eqidx = [26 31 42 43 45 46 87 88 89 ];
M_.mapping.lambda.eqidx = [11 17 20 ];
M_.mapping.qk.eqidx = [11 16 ];
M_.mapping.a.eqidx = [7 33 83 ];
M_.mapping.dn.eqidx = [9 20 34 41 44 80 81 82 83 86 88 89 ];
M_.mapping.sigmaa.eqidx = [33 34 38 ];
M_.mapping.eta.eqidx = [4 27 36 37 ];
M_.mapping.ce.eqidx = [4 27 28 35 40 82 88 ];
M_.mapping.epe.eqidx = [32 ];
M_.mapping.etalag.eqidx = [36 37 ];
M_.mapping.dc.eqidx = [3 80 87 ];
M_.mapping.dcl.eqidx = [2 3 86 ];
M_.mapping.dcc.eqidx = [81 ];
M_.mapping.dce.eqidx = [82 ];
M_.mapping.dcr.eqidx = [87 ];
M_.mapping.dz.eqidx = [83 84 ];
M_.mapping.Edz.eqidx = [84 ];
M_.mapping.dy.eqidx = [88 ];
M_.mapping.di.eqidx = [89 ];
M_.mapping.Phik.eqidx = [9 10 11 ];
M_.mapping.Phik1.eqidx = [11 14 16 ];
M_.mapping.Fk.eqidx = [11 12 ];
M_.mapping.Fl.eqidx = [17 18 ];
M_.mapping.Jk.eqidx = [11 13 ];
M_.mapping.Jl.eqidx = [17 19 ];
M_.mapping.Phip.eqidx = [8 21 39 ];
M_.mapping.Phip1.eqidx = [20 22 ];
M_.mapping.Phip2.eqidx = [20 23 ];
M_.mapping.rf.eqidx = [26 54 55 56 57 94 ];
M_.mapping.d.eqidx = [39 40 41 ];
M_.mapping.dagg.eqidx = [40 44 ];
M_.mapping.dd.eqidx = [41 42 47 48 ];
M_.mapping.ddn.eqidx = [43 ];
M_.mapping.ddr.eqidx = [42 43 49 ];
M_.mapping.ddagg.eqidx = [44 45 50 51 ];
M_.mapping.ddaggn.eqidx = [46 53 92 ];
M_.mapping.ddaggr.eqidx = [45 46 52 90 ];
M_.mapping.pd.eqidx = [47 48 49 ];
M_.mapping.pdagg.eqidx = [50 51 52 53 ];
M_.mapping.rd.eqidx = [48 ];
M_.mapping.rdr.eqidx = [49 54 ];
M_.mapping.rdagg.eqidx = [51 91 ];
M_.mapping.rdaggn.eqidx = [53 93 97 ];
M_.mapping.rdaggr.eqidx = [52 55 96 ];
M_.mapping.rm.eqidx = [54 56 ];
M_.mapping.rmex.eqidx = [56 ];
M_.mapping.yn2.eqidx = [60 61 ];
M_.mapping.yn3.eqidx = [61 62 ];
M_.mapping.yn4.eqidx = [62 63 ];
M_.mapping.yn5.eqidx = [63 64 ];
M_.mapping.yn6.eqidx = [64 65 ];
M_.mapping.yn7.eqidx = [65 66 ];
M_.mapping.yn8.eqidx = [66 67 ];
M_.mapping.yn9.eqidx = [67 68 ];
M_.mapping.yn10.eqidx = [68 69 ];
M_.mapping.yn11.eqidx = [69 70 ];
M_.mapping.yn12.eqidx = [70 71 ];
M_.mapping.yn13.eqidx = [71 72 ];
M_.mapping.yn14.eqidx = [72 73 ];
M_.mapping.yn15.eqidx = [73 74 ];
M_.mapping.yn16.eqidx = [74 75 ];
M_.mapping.yn17.eqidx = [75 76 ];
M_.mapping.yn18.eqidx = [76 77 ];
M_.mapping.yn19.eqidx = [77 78 ];
M_.mapping.yn20.eqidx = [78 79 ];
M_.mapping.slope.eqidx = [79 ];
M_.mapping.fen.eqidx = [59 ];
M_.mapping.cf.eqidx = [90 ];
M_.mapping.dr.eqidx = [91 ];
M_.mapping.cfn.eqidx = [92 ];
M_.mapping.drn.eqidx = [93 ];
M_.mapping.drrf.eqidx = [94 ];
M_.mapping.dry.eqidx = [95 ];
M_.mapping.rmagg.eqidx = [55 57 ];
M_.mapping.rmexagg.eqidx = [57 ];
M_.mapping.Erd.eqidx = [96 ];
M_.mapping.Erdn.eqidx = [97 ];
M_.mapping.en.eqidx = [33 34 37 38 ];
M_.mapping.em.eqidx = [25 ];
M_.mapping.ee.eqidx = [35 ];
M_.mapping.eeta.eqidx = [37 ];
M_.mapping.ed.eqidx = [54 55 ];
M_.mapping.esigma.eqidx = [38 ];
M_.static_and_dynamic_models_differ = false;
M_.has_external_function = false;
M_.state_var = [2 6 7 8 9 10 11 16 17 18 23 24 25 26 27 28 30 49 50 51 58 59 ];
M_.exo_names_orig_ord = [1:6];
M_.maximum_lag = 1;
M_.maximum_lead = 1;
M_.maximum_endo_lag = 1;
M_.maximum_endo_lead = 1;
oo_.steady_state = zeros(97, 1);
M_.maximum_exo_lag = 0;
M_.maximum_exo_lead = 0;
oo_.exo_steady_state = zeros(6, 1);
M_.params = NaN(64, 1);
M_.endo_trends = struct('deflator', cell(97, 1), 'log_deflator', cell(97, 1), 'growth_factor', cell(97, 1), 'log_growth_factor', cell(97, 1));
M_.NNZDerivatives = [366; 488; 436; ];
M_.static_tmp_nbr = [20; 9; 0; 0; ];
close all
M_.params(1) = 0.994;
delta = M_.params(1);
M_.params(2) = 2;
psi = M_.params(2);
M_.params(3) = 10;
gamma = M_.params(3);
M_.params(6) = 1;
lbar = M_.params(6);
M_.params(9) = 0.9;
alphac = M_.params(9);
M_.params(10) = 5;
phic = M_.params(10);
M_.params(4) = 0.5;
xil = M_.params(4);
M_.params(11) = 6;
nu = M_.params(11);
M_.params(12) = 0.33;
alpha = M_.params(12);
M_.params(13) = 0.02;
deltak = M_.params(13);
M_.params(14) = 80;
phip = M_.params(14);
M_.params(15) = 4.8;
xik = M_.params(15);
M_.params(16) = 0.95;
rhon = M_.params(16);
M_.params(24) = 0.983;
rhoa = M_.params(24);
M_.params(25) = (-0.4471);
a0 = M_.params(25);
M_.params(8) = 0.0046;
mu = M_.params(8);
M_.params(23) = 1;
phia = M_.params(23);
M_.params(17) = 0.01836055554715054/sqrt(1-M_.params(24)^2)*sqrt(1-M_.params(16)^2)/0.31224989991992;
phin = M_.params(17);
M_.params(27) = sqrt(1-M_.params(24)^2)*0.09803624925060343;
sigmaa0 = M_.params(27);
M_.params(26) = 0.997;
rhosigma = M_.params(26);
M_.params(29) = 0.005;
sigmatarget = M_.params(29);
M_.params(31) = 1;
corr_vol = M_.params(31);
M_.params(30) = M_.params(31)*M_.params(29)/M_.params(27);
rhosigma_a = M_.params(30);
M_.params(28) = sqrt(max(M_.params(29)^2-M_.params(30)^2*M_.params(27)^2,0));
sigmasigma = M_.params(28);
M_.params(32) = 0.08;
sigmad = M_.params(32);
M_.params(35) = 0.6;
rhoe = M_.params(35);
M_.params(36) = 0.99;
rhoeta1 = M_.params(36);
M_.params(37) = 0;
rhoeta2 = M_.params(37);
M_.params(33) = 0.028*sqrt(1-M_.params(35)^2);
sigmae = M_.params(33);
M_.params(34) = 0.32*M_.params(10)/(M_.params(10)-1)*sqrt(1-M_.params(36)^2*(1+M_.params(37))/(1-M_.params(37))-M_.params(37)^2);
sigmaeta = M_.params(34);
M_.params(38) = 0.02*sign(M_.params(10)-1);
deltaa = M_.params(38);
M_.params(39) = 0.3333333333333333;
ecap = M_.params(39);
M_.params(40) = 0.7;
rhor = M_.params(40);
M_.params(41) = 1.5;
rhopi = M_.params(41);
M_.params(42) = 0.1;
rhoy = M_.params(42);
M_.params(43) = 0.0032;
sigmam = M_.params(43);
M_.params(44) = 0.02955880224154443;
pi_ss = M_.params(44);
M_.params(45) = M_.params(44)-log(M_.params(1)*exp(M_.params(8))^((-1)/M_.params(2)));
yn1_ss = M_.params(45);
M_.params(18) = (exp(M_.params(8))-1+M_.params(13))^(1/M_.params(15));
a1k = M_.params(18);
M_.params(19) = (exp(M_.params(8))-1+M_.params(13))*1/(1-M_.params(15));
a2k = M_.params(19);
rhofx = 1-0.035^2/(2*0.4162^2); 
xifx = sqrt(0.4162^2*(1-rhofx^2));
M_.params(7) = 0.3333333333333333;
l_ss = M_.params(7);
M_.params(50) = ((M_.params(1)^(-1)*exp(M_.params(8)*1/M_.params(2))-(1-M_.params(13)))/(M_.params(12)*(1-1/M_.params(11))*exp((1-M_.params(12))*M_.params(25))))^(1/(M_.params(12)-1))*M_.params(7);
k_ss = M_.params(50);
M_.params(46) = exp((1-M_.params(12))*M_.params(25))*M_.params(50)^M_.params(12)*M_.params(7)^(1-M_.params(12));
y_ss = M_.params(46);
M_.params(51) = (exp(M_.params(8))-1+M_.params(13))*M_.params(50);
i_ss = M_.params(51);
M_.params(47) = M_.params(46)-M_.params(51);
cc_ss = M_.params(47);
M_.params(48) = M_.params(47);
ce_ss = M_.params(48);
M_.params(61) = (1-M_.params(9))/M_.params(9);
pe_ss = M_.params(61);
M_.params(49) = (M_.params(9)*M_.params(47)^((M_.params(10)-1)/M_.params(10))+(1-M_.params(9))*M_.params(48)^((M_.params(10)-1)/M_.params(10)))^(M_.params(10)/(M_.params(10)-1));
c_ss = M_.params(49);
M_.params(52) = M_.params(1)*exp(M_.params(8))^((-1)/M_.params(2));
m_ss = M_.params(52);
M_.params(53) = (-log(M_.params(52)));
r_ss = M_.params(53);
M_.params(54) = (-log(M_.params(52)));
rf_ss = M_.params(54);
M_.params(56) = M_.params(12)*(1-1/M_.params(11))*M_.params(46)/M_.params(50);
Fk_ss = M_.params(56);
M_.params(55) = M_.params(46)*(1-1/M_.params(11))*(1-M_.params(12))/M_.params(7);
Fl_ss = M_.params(55);
M_.params(58) = (-(M_.params(12)/M_.params(11)))/M_.params(50);
Jk_ss = M_.params(58);
M_.params(57) = (-((1-M_.params(12))/M_.params(11)))/M_.params(7);
Jl_ss = M_.params(57);
M_.params(59) = M_.params(55);
w_ss = M_.params(59);
M_.params(64) = (M_.params(9)^M_.params(10)+(1-M_.params(9))^M_.params(10)*M_.params(61)^(1-M_.params(10)))^(1/(1-M_.params(10)));
ph_ss = M_.params(64);
M_.params(5) = M_.params(46)*(1-M_.params(12))*(M_.params(11)-1)/M_.params(11)*((M_.params(6)-M_.params(7))/M_.params(49))^(1/M_.params(4))/(M_.params(7)*M_.params(64));
omega = M_.params(5);
M_.params(62) = M_.params(51)/M_.params(50);
Phik_ss = M_.params(62);
M_.params(63) = 1;
Phik1_ss = M_.params(63);
M_.params(60) = 1;
qk_ss = M_.params(60);
M_.params(20) = 2;
philev = M_.params(20);
M_.params(22) = log(exp(M_.params(8))*M_.params(52)/(1-exp(M_.params(8))*M_.params(52)));
pd_ss = M_.params(22);
M_.params(21) = exp(M_.params(22))/(1+exp(M_.params(22)));
rhopd = M_.params(21);
steady;
%
% SHOCKS instructions
%
M_.exo_det_length = 0;
M_.Sigma_e(1, 1) = 1;
M_.Sigma_e(2, 2) = 1;
M_.Sigma_e(3, 3) = 1;
M_.Sigma_e(4, 4) = 1;
M_.Sigma_e(5, 5) = 1;
M_.Sigma_e(6, 6) = 1;
options_.k_order_solver = true;
options_.irf = 0;
options_.order = 3;
options_.periods = 20000;
options_.pruning = true;
var_list_ = {};
[info, oo_, options_] = stoch_simul(M_, options_, oo_, var_list_);
save('model_file_results.mat', 'oo_', 'M_', 'options_');
if exist('estim_params_', 'var') == 1
  save('model_file_results.mat', 'estim_params_', '-append');
end
if exist('bayestopt_', 'var') == 1
  save('model_file_results.mat', 'bayestopt_', '-append');
end
if exist('dataset_', 'var') == 1
  save('model_file_results.mat', 'dataset_', '-append');
end
if exist('estimation_info', 'var') == 1
  save('model_file_results.mat', 'estimation_info', '-append');
end
if exist('dataset_info', 'var') == 1
  save('model_file_results.mat', 'dataset_info', '-append');
end
if exist('oo_recursive_', 'var') == 1
  save('model_file_results.mat', 'oo_recursive_', '-append');
end


disp(['Total computing time : ' dynsec2hms(toc(tic0)) ]);
if ~isempty(lastwarn)
  disp('Note: warning(s) encountered in MATLAB/Octave code')
end
diary off
